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Topological Characterization of Extended Quantum Ising Models 

G. Zhang and Z. Son^ 

We show that a class of exactly solvable quantum Ising models, inclnding the transverse-field Ising 
model and anisotropic XY model, can be characterized as the loops in a two-dimensional auxiliary 
space. The transverse-field Ising model corresponds to a circle and the XY model corresponds 
to an ellipse, while other models yield cardioid, limacon, hypocycloid, and Lissajous curves etc. 
It is shown that the variation of the ground state energy density, which is a function of the loop, 
experiences a nonanalytical point when the winding number of the corresponding loop changes. The 
winding number can serve as a topological quantum number of the quantum phases in the extended 
quantum Ising model, which sheds some light upon the relation between quantum phase transition 
and the geometrical order parameter characterizing the phase diagram. 

PACS numbers: 03.65.Vf, 75.10.Jm, 05.70.Fh, 02.40.-k 


Introduction. Characterizing the quantum phase tran¬ 
sitions (QPTs) is of central significance to both con¬ 
densed matter physics and quantum information science. 
QPTs occur only at zero temperature due to the com¬ 
petition between different parameters describing the in¬ 
teractions of the system. A quantitative understanding 
of the second-order QPT is that the ground state un¬ 
dergoes qualitative changes when an external parameter 
passes through quantum critical points (QCPs). 

There are two prototypical models, Bose-Hubbard 
model and transverse-field Ising model, based on which 
the concept and characteristic of QPTs can be well 
demonstrated. However, among the two, only the 
transverse-field Ising model is exactly solvable [I| , so as to 
be a unique paradigm for understanding the QPTs. Re¬ 
cently, more attention has been paid to theoretical stud¬ 
ies of exactly solvable quantum spin models involving 
nearest-, next-nearest-neighbor interactions, and multi¬ 
ple spin exchange models, etc 1M3- Those models 
are closer to real quasi-one-dimensional magnets [TTI - [l^ 
comparing to standard ones with only nearest-neighbor 
couplings. Furthermore, it has been shown that quantum 
spin models can be simulated in artificial quantum sys¬ 
tem with controllable parameters. Quantum simulation 
of spin chain can be experimentally realized through neu¬ 
tral atoms stored in an optical lattice [13, [T3| , trapped 
ions [TB - [^ and NMR simulator [2^. This system often 
serves as a test bed for applying new ideas and methods 
to quantum phase transitions. 

A fundamental question is whether QPTs in Ising 
model can have a connection to some topological char¬ 
acterizations. It is interesting to note in this context 
that some simple Ising models have been found to exhibit 
topological characterization [23 - [^ . The purpose of the 
present work is to shed some light upon the relation be¬ 
tween QPTs and a geometrical parameter characterizing 
the phase diagram, through the investigation of a class 
of quantum Ising models. 

In this work, we present an extended quantum Ising 
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model, which includes an additional three-body interac¬ 
tion. It can be exactly solved by the routine procedure, 
taking the Jordan-Wigner and pseudo-spin transforma¬ 
tions. Based on the exact solution, we investigate the 
QPT in this model. We introduce a global order param¬ 
eter, which is the winding number for the loop specifying 
to a set of coupling constants, in an auxiliary space. The 
ground state energy density can be a function of the loop 
and its variation experiences a nonanalytical point when 
the winding number of the corresponding loop changes. 
Then the relation between QPTs and the geometrical or¬ 
der parameter is established. 

Extended Ising model and solutions. We start our anal¬ 
ysis from the one-dimensional Ising model, which has the 
Hamiltonian 
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where (t“, for a = x,y,z, are the usual Pauli matrices, 
and periodic boundary conditions are assumed. Com¬ 
paring with the customary anisotropic XY model, there 
are additional three-site interactions and 

which have the following two implications: 
it can be either regarded as the conditional anisotropic 
XY-type coupling between next-nearest-neighbor spins, 
or conditional action of transverse field. The ground 
state phase diagram and correlation functions for this 
spin model have been studied 40 years ago 0. In the 
case of 5 = 7 = (5 = 0, the correlation function has been 
obtained Q. In addition, other types of Hamiltonians 
which contain three-body interactions were also investi¬ 
gated [30l - l33| . 

We will see that this model can be exactly solvable 
in a simple way by the similar procedure for the sim¬ 
ple transverse-field Ising model [I|, [SJ, [3J|. For the sake 
of simplicity, we only concern the case of even N, the 
conclusion is available for the case of odd N in the ther¬ 
modynamic limit. As the same procedure performed in 
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solving the Hamiltonian without the additional term, we 
take the Jordan-Wigner transformation [l[ 

= ( 2 ) 

<t' = - n (l - 2c,'«) (c, + cj) . (2) 

l<3 

to replace the Pauli operators by the fermionic operators 
Cj. We note that the parity of the number of fermions 
is a conservative quantity and then the Hamiltonian ([T]) 
can be written in the form 

"=(Th-)' W 


where 

H+ =H- -2 
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and 


H- = 
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are corresponding reduced Hamiltonians in the invariant 
subspaces with even and odd number of fermions. Here 
H'^ represents a fermionic ring threaded by a half of the 
flux quantum. In the following, we will focus on H'^ since 
the ground state has even parity for any values of param¬ 
eters. Similarly, can be diagonalized by Fourier and 
pseudo-spin transformations. Taking the Fourier trans¬ 
formation 


where k = 27r (m -I- 1/2) /N, m = 0,1, 2,..., TV — 1, the 
Hamiltonian iJ'*' can be expressed as a compact form 



FIG. 1. (Color online) Schematically illustration of the equiv¬ 
alent Hamiltonian ©, which represents a system of spin en¬ 
semble in a monopole field. In the thermodynamic limit, the 
ground state energy density becomes an integration corre¬ 
sponding to a loop, (b) Schematics of two loops I and I' de¬ 
scribed by the parametric equations involving vectors (k) 
and (fe) -I- 51^ (k), respectively. Here (fei) represents 
an arbitrary point on the I, while (fci) + Sr (fci) repre¬ 
sents the corresponding point on the I'. The red arrow indi¬ 
cates 51^ (ki), while the blue arrow indicates the correspond¬ 
ing unitary vector f (fci). The inner product between them 
51^ (ki) ■ r (ki) contributes to the variation of Sg. The loop 
I passes the origin at point 7^ (fc) with k = ko- The corre¬ 
sponding unitary vector f (feo) is indefinite, which is denoted 
as a solid blue circle. The indefiniteness of 51^ {ko) ■ r (ko) 
witnesses the QPT as well as the topological change of the 
loop: enclosing the origin or not. 


loop of 1^ (fc) defined in Eq. ([5]). In the following ar¬ 
gument, we do not restrict the shape of the loop. The 
obtained result is valid for an arbitrary loop, which is 
schematically illustrated in [IJa). The Hamiltonian ([5]) is 
easy to be diagonalized by aligning all spins with the local 
magnetic field, which is the essential of the Bogoliubov 
transformation. In the thermodynamic limit, the ground 
state energy density can be expressed by an integration 


TT 

— TT 


i7+=4^7^(fc).^fe, (8) 

k>0 

~r (k) = (0, ay sin k + bS sin 2fc, a cos k + b cos 2k — g), (9) 

which represents a set of pseudo spins in a two- 

dimensional magnetic field . The pseudo spin is defined 
as 

■Sfc = (4)^ = CfcC-fe, (10) 

(4cfe + cLfeC-fe - l) , (11) 

satisfying the SU(2) algebra [s^,s^,] = ±(5fcfc'S^,, 

[s/", s/l] = 25kk'S^,. It is clear that the equivalent Hamil¬ 
tonian ([5]) represents a system of spin ensemble in a 
monopole field. These spins locate at the points on the 


which corresponds to a loop tracing with the paramet¬ 
ric equation (fc) = (0, a; (fc), y (fc)). In our case, the 
parametric equation has the form 

J a; (/c) = ay sin fc-b 6(5 sin (2A:) 
yy{k)=acosk + hcos{2k)—g ' ' 

in the auxiliary space (a;,y). Then we can use some 
simple loops to represent the ground state of the Ising- 
like model. It offers many types of graphs correspond¬ 
ing different kinds of Ising model. To demonstrate this 
point, we plot several types of graphs in Fig. O It shows 
that two familiar models, transverse Ising model and 
anisotropic XY model, correspond to two simple graphs, 
circle and ellipse, respectively. Rest models connect to 
more complicated graphs. Furthermore, the ground state 
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(d) a > 2 ti, 7 = 1 , 
5 = 1 


(e) a = 26 ,7 = 1 , 
5 = 1 


(f) a < 26 ,7 = 1 , 
5 = 1 
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(j) a = 0,07 7 ^ 0 , (k) 6 = 0 , 65710 , (I) 6 = 0,7 = 0 

5=0 7=0 



FIG. 2. (Color online) several types of graphs 


of each model is naturally connected to a graph individ¬ 
ually. 

Quantum phase transition. In this section, we inves¬ 
tigate the QPT occurring in the extended Ising model 
and its connection to the geometry of the correspond¬ 
ing loops. To this end, we start with the change of 
the ground state energy induced by varying the parame¬ 
ters. In general the QPT driven by the parameters {a} 
(a = a, b, 7 , A, g) can be characterized by the derivative of 
the ground state density with respect to a, degjda which 
experiences a nonanalytic point at the critical point, lead¬ 
ing to the divergence of second derivative of ground state 
energy density. We investigate the signature of the QPT 
in an alternative way: Eg depends on the path of the in¬ 
tegral, being a function of the functions x {k) and y (k). 
The parameters {a} drive the QPT through the change 
of the functions x{k) and y{k), or the loop. In other 
words, one can consider the variations of functions x {k) 
and y (k) instead of the change of the parameters {a}. 
The first variation of the function eg[x,y] is 

TT 

— TT 

(14) 

where f (fc) = A/ |A| is the unitary vector of A (fc). It 
indicates that the variation Ssg is the summation of the 
path shifts Sl^ (fc) along the direction of f{k). We are 


interested in the case of the loop crossing the origin. At 
the origin, the unitary vector r (k) is indefinite, which 
leads to an indefinite contribution to the variation Ssg, 
indicating a nonsmooth point. It is a signature of the 
QPT associated with a topological change in the loop of 
the integration. So far, we do not specify the shape of 
the loop and how the loop is deformed. In our case, the 
variation (5 A (fc) arises from the continuous change of the 
parameters {a}. Then we have 

jA (fc) = y ^^-^da (15) 

da 

a 


or explicitly 

f 5x (fc) = sin k (yda -I- ady) -I- sin {2k) {bdX + Xdb) 

\ Sy (k) = cos kda + cos {2k) db — dg 

(16) 

Considering the case with a = 7 = l ,6 = A = 0as 
an example, the Hamiltonian o reduces to the simplest 
transverse field Ising model 

N 

ifising = y (ayV+ga^^), (17) 

i=i 

the ground state energy density of which corresponds to 
a circle of the equation 

x^ + {y + 9)^ = ^- (18) 

The variation Ssg from the case with 5 = ±1 is readily 
expressed as 


TT 

6eg = -^ J [f{k)-j]dk, (19) 

— TT 

where j denotes the unit vector of y axis. It results 
in Ssg = dsg/dgdg, which shows that Ssg accords with 
Osg/dg as a witness of QPT. 

We would like to point out that each loop contains 
two characters, geometry (shape and position) and curve 
orientation, which are determined by the corresponding 
parameter equation. To characterize these two features, 
we use a topological quantity, winding number, which is 
a fundamental concept in geometric topology and widely 
used in various areas of physics |35l438l |. The winding 
number of a closed curve in the auxiliary xy-plane around 
the origin is defined as 

■A = ^ ^ ^ (ydx - xdy ), (20) 

which is an integer representing the total number of times 
that the curve travels clockwise around the origin. Then 
we establish the connection between the QPT and the 
switch of the topological quantity. 

Topological quantum number. We calculate the wind¬ 
ing numbers for various typical cases corresponding to 
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FIG. 3. (Color online) The winding numbers for various typ¬ 
ical cases corresponding graphs in an auxiliary space with 
different topologies. Label f (4-) denotes the increase (de¬ 
crease) of the parameters, which induce the transition be¬ 
tween graphs. 


graphs in an auxiliary space with different topologies. In 
Fig. [31 we plot the graphs of the ground state for typical 
cases. The corresponding winding number and the rela¬ 
tions between graphs are presented. In each group, the 
first graph is clockwise. The examples show that there 
are five possible winding numbers ±2, ±1, and 0, which 


represent five different phases. 

To demonstrate the characteristics of these phases, 
we consider five typical cases, which correspond to the 
ground states of systems with parameters in fol¬ 
lowing limits: (i) & —>■ oo and d = — 1 , the reduced 
Hamiltonian is /i _2 = (ii) a —>■ 

oo and 7 = -1 and /i_i = (iii) 

g —>■ 00 , ho = (iv) a —>■ 00 and 7 = 

1, hi = (v) 6 -)> 00 and (5 = 1, 

/i 2 = N~^ corresponding ground 

states IGo.ibi) of even-number flip subspace, represented 
in position space, are readily obtained as 

|Go) = n 1 ^).- ’ ( 21 ) 

3 

iG.)=^(n I/-),' n i/>i+n i/). n <^2) 

jGe j^o j^e j^o 

iG-.>=^(n i\), n IN),+n in), n in>,). m 

j^e j^o j£e j^o 

where cr| li)^-(|t)j) = - 14-)^ (lt)j), cr“ |/’)^-(|^/)^-) = 

IA (-|^).)’ ^?I\),(N),) = l\),(-|\),)( ande 

and o denote the even and odd number of sites, respec¬ 
tively. The ground states of h ±2 are obtained from Eq. 
(I5|) and expressed in an auxiliary space as 

\G'± 2 ) = n (±isinfc|t)fe + cosfc|4.)^), (24) 

k>0 

where |t)fc and 14.);, are eigenstates of pseudo-spin oper¬ 
ator si with 2sl |t)fe (|4.)fc) = lt)fe (- I4.)fe)- It is a little 
complicated to express states |G± 2 ) in the position space 
in a simple form. Here we only give the expression for 
N = An{nGN) [H, 


|G± 2 ) = 2 -i^- 2)/2 

n 


N/2 

E(±i) 

3=0 




^ ^ ^ \ ^ni 


(25) 


where 1^) = n;ii It), is the saturate ferromagnetic state. 
As an example, the ground states for 4-site systems h ±2 
are explicitly 

|G^2="> = l/2(|t)i lt)2 lt)3 lt)4 T |t)i lt)2 lt)3 lt)4 

T|t)ilt)2lt)3lt)4-|t)ilt)2lt)3lt)4)-(26) 


We employ the expected value of operators hp, 
{G\\hp\G\), as local order parameters to characterize 
the ground states IGA {p,X = ±2, ±1,0). By using the 
similar analysis in [39j, we have 


{Gx\hp\Gx) = -6xp. (27) 


It indicates that the five ground states |Ga) are in five dif¬ 
ferent phases. Then the winding number can be reliable 
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topological quantum number to distinguish the quantum 
phases. 

Conclusion. In summary, a class of exactly solvable 
quantum Ising models presented in this paper have ob¬ 
vious topological characterization and indicate the ex¬ 
istence of a topological quantum number, which is the 
winding number for the loop in a two-dimensional auxil¬ 
iary space and describes the quantum phases in the ex¬ 
tended quantum Ising model. This finding reveals the 
connection between QPT and the geometrical order pa¬ 


rameter characterizing the phase diagram for a more gen¬ 
eralized spin model, which will motivate further investi¬ 
gation. 
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